Interacting generalized ghost dark energy in non-flat universe 
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We investigate the generalized QCD ghost model of dark energy in the framework of Einstein 
gravity. First, we study the non-interacting generalized ghost dark energy in a flat FRW background. 
We obtain the equation of state parameter, wd = p/p, the deceleration parameter, and the evolution 
equation of the generalized ghost dark energy. We find that, in this case, wd cannot cross the 
phantom line (wd > —1) and eventually the universe approaches a de-Sitter phase of expansion 
(wd — ► —1). Then, we extend the study to the interacting ghost dark energy in both flat and non- 
flat FRW universe. We find out that the equation of state parameter of the interacting generalized 
ghost dark energy can cross the phantom line (wd < —1) provided the parameters of the model are 
(N. chosen suitably. We also discuss the values of the deceleration parameter in all cases and present a 

numerical discussion of this parameter in different phases of the universe expansion. 
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I. INTRODUCTION 



A complementary cosmological observations confirm that our universe is currently undergoing a phase of accelerated 
expansion [1]. Within the framework of general relativity (GR), the acceleration can be explained by assuming a bizarre 
type of energy with a negative pressure called "dark energy" (DE). A plenty of models have been proposed to explain 
the acceleration of the universe expansion. One can refer to [HIH for a review of DE models. On the other side, many 
people believe in a modification of gravity, seeking an explanation for the late time acceleration. According to this 
Ch . idea the acceleration will be a part of the universe expansion and does not need to invoke any kind DE component. 
As examples of this approach one can look at Refs. (3-0]. 

DE models are mainly based on introducing a new degree of freedom or a new parameter which provide the needed 
conditions for explanation of the cosmic acceleration (see e.g. Q and references therein). However, it would be nice 
^ ' to resolve the DE puzzle without presenting any new degree of freedom or any new parameter in the theory. One 
of the successful and beautiful theories of modern physics is Quantum Chromodynamics (QCD) which describes the 
strong interaction in nature. However, resolution of one of its mysteries, the U(l) problem, has remained somewhat 
unsatisfying. Veneziano ghost field explained the U(l) problem in QCD [9]. Vacuum energy of the ghost field can be 
used to explain the time varying cosmological constant in a spacetime with nontrivial topology, since the ghost filed 
has no contribution to the vacuum energy in the Minkowskian spacetime The ghost vacuum energy density is 
proportional to Aq CD H, where Aqcd is QCD mass scale and H is the Hubble parameter [ll|. It is well-known that 
the cosmological constant model of DE suffers the coincidence and the fine tuning problems. However, with correct 
• • ■ choice of Aqcd, the ghost dark energy (GDE) model does not encounter the fine tuning problem anymore [lOL 

Phenomenological implications of GDE model were discussed in [12]. In [r| the GDE in a non-fiat universe in the 
presence of interaction between DE and dark matter was explored. The instability of GDE model against perturbations 
was studied in [14j . It was argued that the perfect fluid for GDE is classically unstable against perturbations. Other 
features of GDE model have been investigated in Refs. p^l - [r7| . 

In all above references f [Iol - l2l| ) the GDE was assumed to have the energy density of the form pr> = aH, while, in 
general, the vacuum energy of the Veneziano ghost field in QCD is of the form H + 0(H 2 ) (22|. This indicates that 
in the previous works on the GDE model, only the leading term H has been considered. Motivated by the argument 
given in 23[, one may expect that the subleading term H 2 in the GDE model might play a crucial role in the early 



evolution of the universe, acting as the early DE. It was shown |24| that taking the second term into account can give 
better agreement with observational data compared to the usual GDE. Hereafter we call this model as the generalized 
ghost dark energy (GGDE) and our main task in this paper is to investigate the main properties of this model. In 
this model the energy density is written in the form po = olH + f3H 2 , where fj is a constant. 

Beside the DE component in the universe there exist another unknown component called "dark matter" (DM). Since 
the nature of these two dark components are still mystery and they seem to have different gravitational behavior, 
people usually consider them separately and take their evolution independent of each other. However, there exist 
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several observational evidences which detect signatures of interaction between the two dark components. For example, 
observational evidences provided by the galaxy cluster Abell A586 support the interaction between DE and DM [25| ■ 
Also, it has been shown that interaction can be compatible with SNIa and CMB data [IB]- Considering the interaction 
can also be understood from theoretical point of view. When we assume that the DM and the DE components evolve 
separately it means that both of them are conserved independently. Any conservation results from a symmetry in the 
Lagrangian level and in the absence of such a symmetry it seems acceptable to take an interaction between DE and 
DM. 

Based on the cosmological principle the universe has three distinct geometry, namely open, flat and closed geometry 
corresponding to k = —1, 0, +1, respectively. For a l ong time it was a general belief that universe has a flat (k = 0) 
geometry which mainly based on the inflation theory [2]| ■ With the development of the observational technics people 
found deviations from the flat geometry (28|. For example, CMB experiments [2t| leave the chance for a spatially closed 
universe. In addition to CMB, recently the spatial geometry of the universe was probed by supernova measurements 
of the cubic correction to the luminosity distance |3Q|, where a closed universe is also marginally favored. Also the 
information resulted from WMAP data provides further confidence, showing that a closed universe with positively 
curved space is marginally preferred l_3l]. 

All above reasons, indicate that although people believe in a flat geometry for the universe, but astronomical 
observations leave enough room for considering a nonflat geometry. Also about the interaction between DM and DE 
there are several signals from nature which guides us to let the models explain such a behavior. Based on these 
motivations we would like here to extend the studies on GGDE, to a nonflat FRW spacetime in the presence of an 
interaction term. 

This paper is organized as follows. In the next section, we review the GGDE model in a flat universe. Section HTI1 
is devoted to the interacting GGDE in a flat geometry. In section ITVl we generalize the study to the universe with 
spacial curvature in the presence of interaction between DM and DE. We summarize our results in section IVl 



II. GGDE MODEL IN FLAT UNIVERSE 

We start with studying the non interacting GGDE model in a flat universe. Although this model has been studied 
in [24j |. but we would like to consider the GGDE in a different manner by direct look at the important parameters 
such as the EoS parameter, wd, and the deceleration parameter, q. 

For flat Friedmann-Robertson- Walker (FRW) universe, the first Friedmann equation takes the form 

H 2 = —(p m + PD ), (1) 

where p m and po are, respectively, the energy densities of pressureless DM and DE. The generalized ghost energy 
density is defined as [24| 

PD = aH + l3H 2 , (2) 

where a is a constant of order Aq CD and Aqcd is QCD mass scale, and j3 is also a constant. In the original GDE 
(/3 = 0) with Aqcd ~ WOMeV and H ~ 10 _33 eV , A QCD H gives the right order of magnitude - (3 x !CT 3 cV) 4 for 



the observed DE density In the GGDE, j3 is a free parameter and can be adjusted for better agreement with 
observations. 

We define the dimensionless energy density parameters as 

p m 8irGp m _ Pd _ 87rG(a + (3H) 

Per dH £ p cr 3H 

where p cr — 3H 2 / (8irG). According to these definitions, the first Friedmann equation can be rewritten in terms of 
the energy density parameters as 

n m + n D = i. (4) 

Through this section we consider GGDE in the absence of interaction term, thus the conservation equations read 

p m + 3Hp m = 0, (5) 
p D + 3Hp D (l + w D ) = 0, (6) 



which imply that DM and DE are separately conserved. 
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Taking time derivative of relation ((2J, we obtain 

Pd = H(a + 2/3H). (7) 

Also it can be easily seen that differentiating Eq. (TTJ) with respect to time leads to 

H = -AnGp D (l + u + w D ), (8) 

where u = p m / Pd is the energy densities ratio. Replacing this relation in Eq. (JT)) and also using the continuity 
equation (JSJ), we get 

(1 + w D )[3H - 4nG(a + 2/3H)} = AnG{a + 2(3H). (9) 
Solving the above equation for wd and noticing that u — f2 m /fii) and 

— ( a + 2pH) = -?- + ^, (10) 

we obtain 

WD = n D (2-n D -o> 

where £ = 87r ^ . It is clear that this relation reduces to its respective one in the GDE when £ = [13| . 

In Fig. la we plotted the evolution of u>d versus f2_D. It is easy to see that at the late time where f2_o — > 1, we have 
wd — > — 1) which implies that the GGDE model mimics a cosmological constant behavior. One should notice that 
this behavior is the same as for the original GDE model. This is expected since the subleading term H 2 in the late 
time can be ignored due to the smallness of H and the difference between these two models appears only at the early 
epoches of the universe. From figure (la) we see that wd of the GGDE model cannot cross the phantom divide and 
the universe has a de Sitter phase at the late time. It is important to note that the universe is filled with two dark 
components namely, DM and GGDE. Thus to discuss the acceleration of the universe we should define the effective 
EoS parameter, w c s, as 

Pt Pd rio s 
Weft = — = ; , (12) 

Pt PD + Pm 

where pt and pt are, respectively, the total energy density and the total pressure of the universe. As usual, we have 
assumed the DM is in the form of pressureless fluid (p m =0). Using relation ^ for the spatially flat universe, one can 
find 

w cS = Q d Wd = » - ni ^ D (13) 

I — Md — 4 

Let us now turn to the deceleration parameter which is defined as 

ad H 

* = -^ = -l-^ (14) 
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where a is the scale factor. Using Eq. ([8]) and definition J7d in ([3]) we obtain 



H 3 

— = --n D (i + u + w D ). (is) 



Replacing this relation into (fl4)) . and using (fTTj) we find 

1 3 C-^u 



y ~2 2(^ + n D -2)- (16) 

One can easily check that the deceleration parameter in GDE is retrieved for £ = [13|. We can also take a look at 
the early and the late time behavior for the deceleration parameter of this model. At the early time where Sl^ — ^ 
the deceleration parameter becomes 

which indicates that for £ < 2 the universe is at the deceleration phase at early times while for £ > 2, the universe 
could experience an acceleration phase, the former is consistent with the definition £ = 8 *GP , On the other side, we 
find that at the late time where the DE dominates (fi^ — » 1), independent of the value of the £, we have q = — 1. We 
have plotted the behavior of q in Fig. [Tb. Besides, taking Sl D0 — 0.72 and adjusting £ = 0.01 we obtain q n ~ —0.34 for 
the present value of the deceleration parameter which is in agreement with recent observational data [32j |. Choosing 
the same set of parameters leads to wdo ~ —0.78 and w e so ~ —0.56. One should note that as we already mentioned 
about Wb-i the square term in the GGDE density has a negative contribution in the role of the DE in the universe. 
We mean by negative contribution that arises by taking the square term into account, the evolution of the universe 
will be slowed. For example, the universe will enter the acceleration phase later than the original GDE. This behavior 
is clearly seen in both parts of Fig. [TJ 

At the end of this section we present the evolution equation of the DE density parameter Vlo- To this goal we take 
time derivative of Eq. ([3]), after using relation Qjj = H Ift 5 - as well as Eq. (fT4")l we reach 

xi D (i-n D )vj D . (is) 



d In a 



d\na 

Using Eq. ((TTj) we get 

dn D = (i-n D )(z-n D ) 

d\na 2-Q D -Z ' [ ' 

Once again for the limiting case £ = 0, the above relation reduces to its respective evolution equation for the original 
GDE presented in QJ]. 



III. INTERACTING GGDE IN FLAT UNIVERSE 



In the previous section, the evolution of DE and DM components were discussed separately. Here we would like 
to extend the study to the interacting case, seeking new features of GGDE. In the first look investigating interacting 
models of DE are valuable from two perspective. The first is the theoretical one, which states that we have no 
reason against interaction between DE and DM components. For example, in the unified models of field theory DM 
and DE can be explained by a single scalar field, thus they will be allowed to interact minimally. Besides, one can 
get rid of the coincidence problem by taking into account the interaction term between DM and DE. One can refer 
to [13-113] for detailed discussion. The other feature which motivates us to consider interacting models of DE and DM 
comes from observations. For instance, observational evidences provide by the galaxy cluster Abell A586 support the 



interaction between DE and DM [25J. Thus, there exist enough motivations to consider the GGDE in the presence of 



an interaction term. To this end, we start with the energy balance equations for DE and DM, namely 

p m +3Hp m = Q, (20) 
PD +3H PD {l + w D ) = -Q, (21) 

where Q > represents the interaction term which allows the transition of energy from DE to DM. The form of Q is 
a matter of choice and can be taken as [131 ] 

Q = Sb 2 H( Pm + PD ) = 3b 2 Hp D (l + u), (22) 
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with 6 2 being a coupling constant. Inserting Eqs. (0 and ([22)) in Eq. (|2ip and taking into account u = ^ 2a , we find 



Wd = — 



2 - - £ 



1 



262 



(23) 



At first look one can find that setting 6 — 0, wd reduces to the respective relation in the absence of interaction 
obtained in Eq. (fTTj) . When £ = the result recovers those in [13] for original GDE. The first interesting point 
about the EoS parameter of the GGDE is that in the interacting case independent of the interaction parameter, 6, for 
< £ < 1, wd can cross the phantom line in the future where Qd — ► 1. At the present time, by choosing £ = 0.03, 
b = 0.15 and ^Ijjq — 0.72, we find that woo = —0.82 and w c so — —0.59 which the latter favored by observations. 
One can easily check that for a same coupling constant these values for the original GDE are wdo = —0.83 and 
Weffo = —0.60 which clearly show that the square term in the energy density of the GGDE slow down the evolution 
of the universe compared to the original GDE model. For a better insight we have plotted wd against O^ in Fig. [5^- 
This value for coupling constant, 6, in the figure is consistent with recent observations (38[. It is worth mentioning 
that at the late time where — !> 1 the effective EoS parameter approaches less than —1, i.e. w c g < —1, which 
reminds a super acceleration for the universe in the future. Next we take a look at the deceleration parameter in the 
presence of an interaction term. Substituting (fT5)l in (fT4")l and using (|2"3"|) yields 

262 



n 



D 



1 



±_ 



(24) 



* 2 2 (2 - n D - o 

Once again it is clear that setting 6 = 0, the respective relation in the previous section is retrieved. When £ = 
the result of [l3[ is recovered. For the set of parameters (£ = 0.03,6 = 0.15, SIdo = 0.72), we find that according to 
the GGDE the universe enters the acceleration phase at = 0.48 while this transition happens earlier for the GDE 
model. This point is clear from FigHJj. The present value of the deceleration parameter for the interacting GGDE 
model is go = —0.38 which is consistent with observations (32| . 

Finally, we would like to obtain the evolution equation of DE in the presence of interaction. First we take the time 
derivative of ^ and obtain 



n = n 



P H 



Using relation (|2"Tj) as well as (|T5|) , it is a matter of calculation to show 



din a 



= 3fi; 



1 - 



D 



2-0, 



26 2 




b 2 









(25) 



(26) 



In the limiting case £ = the equation of motion of interacting GDE is recovered [i± 



IV. INTERACTING GGDE IN NON-FLAT UNIVERSE 



Flatness problem in standard cosmology was resolved by considering an inflation phase in the evolution history of 
the universe. Following this theory it became a general belief that our universe is spatially flat. However, later it was 




shown that exact flatness is not a necessary consequence of inflation if the number of e- foldings is not very large [39j . 
So it is still possible that there exist a contribution to the Friedmann equation from the spatial curvature, though 
much smaller than other energy components according to observations. Thus, theoretically the possibility of a curved 
FRW background is not rejected. In addition, recent observations support the possibility of a non-flat universe and 
detect a small deviation from k = |40j|. For example evidences from CMB and supernova measurements of the 
cubic correction to the luminosity distance favor a positively curved universe [U |42[. Besides, some exact analysis 
of the WMAP data support a closed geometry for our universe [43j]. Furthermore, the parameter fi^ represents the 
contribution to the total energy density from the spatial curvature and it is constrained as —0.0175 < £lk < 0.0085 
with 95% confidence level by current observations [44]. Our aim in this section is to study the dynamic evolution of 
the GGDE in a universe with spatial curvature. As we discussed in the introduction a closed universe is marginally 
favored. We start with the first Friedmann equation in a non-flat universe which is written as 

H 2 + ^ = ^-A Pm + Pn), (27) 

where k is the curvature parameter with k = —1,0, 1 corresponding to open, flat, and closed universes, respectively. 
Taking the energy density parameters ([3]) into account and defining the energy density parameter for the curvature 
term as flk — k/(a 2 H 2 ), the Friedmann equation can be rewritten in the following form 

i + n k = n m + n D . (28) 

Using the above equation the energy density ratio becomes 

Pm 1 + — /r .„ N 

u = — = o~ = o • ( 29 ) 

PD ' 'D "D 



The second Friedmann equation reads 



H = -4nG(p + p) + ^, (30) 



while the time derivative of GGDE density is 

Pd = H(a + 2/3H). (31) 



Inserting Eq. ([30]) into (|3Tj) and combining resulting relation with conservation equation for DE component (|21j). after 
using (|22p and (j2"9"]l , we find the EoS parameter of interacting GGDE in non-flat universe 



From the second Friedmann equation, (|30l) . one can easily obtain 
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H 3 

— = - fl k + -Q. D [1 + u + w D ], 
and therefore the deceleration parameter in a non-fiat background is obtained as 



(33) 



H 3 

i - -=5 = -i - n* + -n D [i + u + w D ]. 



Substituting Eqs. (gSJ) and ([32) in (EU we obtain 
^I-lO^ 3» g 

q =2 {1 + nk) - 2(2-n D -o 



(34) 



(35) 



In a non-flat FRW universe, the equation of motion of interacting GGDE is obtained following the method of the 
previous section. The result is 



D 



din a 



30; 



~3~ 



1 - Q, 



fi D 



-(l + fife) 



(36) 



In the limiting case flk = 0, the results of this section restore their respective equations in flat FRW universe derived 
in the previous sections, while for £ = the respective relations in [13[ are retrieved. The evolutions of wd and q 
against fljj for a non-flat interacting GGDE and GDE models are plotted in FigJ3J Let us explore different features 
of GGDE in non-flat universe by a numerical study. First of all we study the EoS parameter of the GGDE in the 
future where flo — > 1- In this case, taking £ = 0.1, b = 0.15 and Q k = 0.01 leads to wd = —1.05 which indicates that 
the GGDE is capable to cross the phantom line in the future. The present stage of the universe can be achieved by 
the same set of parameters but flD = 0.72. In such a case we see that wdq = —0.78 while the effective EoS parameter 
becomes w e go = —0.6 which is consistent with observations. The deceleration parameter of the model can also be 
obtained which is in agreement with observational evidences. For example, for the above choice of parameters one 
finds qo = —0.34 [32| . Transition from deceleration to the acceleration phase, in the interacting non-flat case, take 
place at fljj — 0.52. 



V. SUMMARY AND DISCUSSION 



Veneziano ghost model was presented to solve U(l) problem in QCD theory. People usually prefer to handle the 
DE problem in physics using degree's of freedom already exist. GDE is a prototype of these models which discusses 
the acceleration of the universe and is based on the mentioned Veneziano ghost field. This model of DE can address 
the fine tuning problem [l3{. An extended version of this model called GGDE was recently proposed by Cai et. al., 
[24j |. seeking a better agreement with observations. 

In this paper we explored some features of GGDE in both flat and non-flat FRW universe in the presence of an 
interaction term between the two dark components of the universe. In section II, we discussed the GGDE in a flat 
FRW background. We found that the EoS parameter approaches —1 which is the same as the cosmological constant. 
The next section was devoted to the interacting GGDE in flat geometry. An interesting feature which we found was 
the capability of crossing the phantom line in this case. This behavior is also seen in the last section for interacting 
GGDE in a universe with spatial curvature. 

We also investigated the deceleration parameter in different cases. For flat universe in the presence of interaction 
between DE and DM we found that with fl^ = 0.72 for present time and choosing £ = 0.01 and b — 0.15, we reach 
q = —0.38 which shows a consistent value in comparison with observations (32l |. Also this parameter has a consistent 
behavior in all cases of interest. 

In summary the main result found through this paper is that in the GGDE model, there is a delay in different 
epoches of the cosmic evolution in comparison with original GDE model. For instance we saw that the deceleration 
parameter for the present time in the flat noninteracting GDE was —0.34 while for GGDE, the same parameter was 
obtained —0.33 for £ = 0.03. This point is completely seen from the figures presented in the paper. This result were 
also pointed in [24| due to the negative contribution of the square term in the energy density of GGDE. 
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